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Abstract 



Systems of spin 1, such as triplet pairs of spin-1/2 fermions (like ortho- 
hydrogen nuclei) make useful three-terminal elements for quantum com- 
putation, and when interconnected by qubit equality relations are univer- 
^ ' sal for quantum computation. This is an instance of quantum-statistical 

I computation: some of the logical relations of the problem are satisfied 

. identically in virtue of quantum statistics, which takes no time. We show 

heuristically that quantum-statistical ground-mode computation is sub- 
stantially faster than pure ground-mode computation when the ground 
mode is reached by annealing. 

1 Introduction 

1.1 Premise 

' To-day quantum computation is mostly algorithmical: a stored program controls 

^ , a time-varying Hamiltonian that sends a set of qubits - two-valued quantum 

variables - through a sequence of unitary transformations effecting steps in a 
computation. The quantum speed-up depends essentially on maintaining co- 
herent parallel computation paths during the entire computation, until the final 
measurement. 

' The well-known fragility of such a coherence may make algorithmic quantum 

computation impractical. 

Quantum ground-mode computation could be an interesting alternative. One 
defines a time-independent quantum Hamiltonian proportional to the amount 
of "frustration" (violation) of the logical relations of the problem to be solved. 
Then any ground mode solves the problem. 

An incoherent superposition of ground modes also solves the problem. This 
should greatly reduce the problem of decoherence. 

In quantum annealing computation, one form of quantum ground-mode com- 
puytation, the quantum computer is brought to minimum energy by a suitable 
coupling with a heat bath of programmed temperature. In quantum adiabatic 
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computation (Farhi ct al.), another form, the computer is prepared in a mini- 
mum of a simpUfied Hamiltonian, which is then adiabatically deformed into the 
problem Hamiltonian, leaving the network at minimum energy. 

To be sure, quantum ground-mode computation is suspected of being math- 
ematically and physically intractable. At any rate, its computer simulations are 
lengthy and seriously limit research on this method; and since the numbers of 
energy levels, energy traps (for annealing computation), and energy- level cross- 
ings or proximities (for adiabatic computation), can grow exponentially with 
problem size, ground-mode computation alone may not yield any speed-up. 

We ameliorate this problem here. We implement part of the Boolean re- 
lations as quantum symmetries expressing particle indistinguishability, using a 
quantum-statistical three-terminal network element. Satisfying such relations 
does not take time. There is a qiiantum speed-up for annealing computation in 
the sense that the speed of relaxation of a hard problem becomes comparable to 
that of an easy problem, which algorithmically can be solved in polynomial time 
by successively eliminating variables. Implementing part of the Boolean rela- 
tions through quantum symmetries due to particle indistinguishability should 
give a speed-up over the case where all relations are implemented through an 
energy function. 

1.2 Outline 

We express the problem to be solved in terms of the computationally-universal 
relations Qx+qy + Qz = 1, called the (sum-1) triode, and = Qy, called the wire. 
Here the q are Boolean variables and -t- denotes arithmetical sum. Quantum- 
statistical computation implements the wires through Hamiltonian terms, and 
the triodes through quantum statistics (see also Castagnoli 1998, Castagnoli et 
al. 1998, CastagnoH & Monti 1999). 

We associate each wire qm — qn with a Hamiltonian term (g™ — qn)"^ whose 
two-fold degenerate ground modes satisfy the wire. We associate each triode 
with a spin pair that satisfies the relation identically in virtue of particle statis- 
tics, without time development. The associated Hamiltonian is 0. Satisfying 
relations implemented by statistics does not take time or energy. 

We show this here for quantum-statistical annealing (qusa) computation, by 
means of a special representation of the computation process. 

The effective Hamiltonian H = + Hr{t) consists of a Hamiltonian iJ^, 
for the wires of the network, and a time-dependent heat-bath coupling Hj. to 
relax the network to its zero point. There is no triode Hamiltonian. 

To estimate the solution time, we introduce a certain non-symmetrized com- 
parison Hamiltonian H' that also includes a coupling to a heat bath and yields 
the actual network Hamiltonian H when statistics - symmetrization - is im- 
posed. H' describes the relaxation process of an easy network, where all the 
Boolean constraints implemented by statistics are removed. The easy network 
may be realized in principle by replacing each pair of identical-indistinguishable 
fermions forming a triode by a pair of distinguishable (e.g. non-identical) 
fermions. But when we continuously project the process governed by H' on 
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the Hilbcrt space satisfying the statistics, we recover the relaxation process of 
the actual network. 

In this way we show that the solution time for the actual network is compa- 
rable with that of the easy network. 

Although the interplay between relaxation and statistics proposed is a well- 
defined physical effect, the computation model based on it still dwells in the 
same abstract conceptual realm as other current literature on quantum ground- 
mode computation (see for example Farhi et al., among others). Our model 
demonstrates that a new form of quantum speed-up is possible in principle, 
leaving implementation problems for the future. 



2 The Boolean problem 

By a triode with qubits QxTQy, Qz we mean the relation Qx + % + = 1; where 
+ denotes arithmetical sum. This is the negation of the FOR (partial OR) 
relation used by Boole and algebraicized by C. S. Peirce [Finkelstein 1996]. It is 
a partial NOR (or partial Scheffer stroke) relation and can be written as Qz = Qx 
PNOR qy . 

We consider a Boolean triode network, one consisting only of T triodes, 

Q = 3T qubits, and W wires between qubits. 
The triode is satisfied by just the values 
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qz 
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(1) 



No value of q^ is defined for qy = = 1\ and similarly for qy and g^. This 

triode is therefore not a gate, since it does not define an input-output function 
relating its variables. We may call this triode a partial gate, since it defines a 
partial function; in fact, three. 
A toy triode network is 

91 — (fi 

Xl — 92 — 92 — IT 



with six qubits qx, q2, qs, q[,q'2, q'3, two triodes gi -|- 92 + 93 = 1 , 9l + 92 + 93 = 1 
and three wires qi = q'l , 92 = 92 : qs = q's- 

The problem is to solve the network assuming that there is at least one 
solution. This problem has practically the same difficulty as checking whether 
the network is satisfiable. It is the "exact cover" problem addressed in quantum 
adiabatic computation (Farhi et al. 2001), and is NP-complete. 
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3 The network model 



3.1 The qubit 

The qubits used in quantum computation differ from the bits of classical com- 
putation in that a qubit q has variables that do not commute with q, and modes 
that are quantum superpositions of the q = and q = 1 modes. Qusa compu- 
tation exploits the quantum nature of the qubit still further. 

3.2 The sum-1 triode 

We model the triode with a triplet proton pair identically fulfilling the triode 
relation q^ + qy + qz = 1 in virtue of statistics, as follows. 

A quantum spin 1/2 provides three anticommuting two- valued variables 
ax, (Ty. (T,r, each taking two values ±1, and subject to the relation ax^yaz = i- 

A hydrogen molecule has a triplet of ground modes and an excited singlet 
mode. Let \a\, \ai be the respective spin vectors of protons 1 and 2 (in units 
with fi = 1). Then s = \((y\ -Va-i) represents the total spin angular momentum. 

Two independent spin vectors a\ and a-i define a total spin s = (a\ + cr2)/2 
with Sz = ±1, 0, and three commuting binary variables qx = 1 — s^, Qy = 1 — Sy, 
qz = 1 ~ si subject to the relation qx + qy + Qz = ^ mod 2. These spins thus 
constitute a gate - implementing the weaker relation qx + Qy + Qz = ^ mod 2 - 
among the three qubits: 
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1 


1 



(3) 



Define the spin quantum number S by 

s^ = S{S + l) (4) 

as usual. The fourth line in this table is the singlet mode S = 0. This gate is the 
(non-universal) complemented-XOR gate. Since each variable is 1 (true) when 
and only when the other two are equal, it may also be called the EQUALS (or 
EQU) gate. 

EQU is a functional relation. It may be solved for any of its qubits, say qx- 

qx = qy + qz + i mod 2. (5) 

Networks of triodes and wires are universal for Boolean computation; it is 
easy to construct NOT from two triodes and NOR from three triodes, with 
several wires. 

Networks of EQU gates and wires are not universal. A Boolean network 
made of just EQU gates and wires is a system of modulo-2 arithmetical equations 
that is quickly solvable. 
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In particular there is always the solution where all qubits are 1 and all triodes 
are singlet. 

When = the total spin is "up-or-down," which is doubly degenerate. 
When = 1 it is "not up-or-down," which is also doubly degenerate. The 
qubits Qx and qy have similar meanings relative to the x and y axes. 

In the singlet mode, Qx = Qy = Qz = ^■ 

If we restrict the system to its triplet mode S = 1, the three commuting 
qubits (qx^Qy^Qz) obey the PNOR relation 

qx + qy + qz = '^ (6) 

by (4). Then "not up-or-down" becomes (approximately) "sideways" and is 
non-degenerate. "Up-or-down" remains doubly degenerate. 

In the comparison network constructed later we drop the restriction to the 
triode relation. All triodes are replaced by EQU gates of truth table (3). 

3.3 The error metric 

We define an "error metric" measuring the distance of the network from a 
solution. Prom now on we use the following notation. We index the network 
wires with ui = 1,. . . ,W and the triodes with r = 1, ...,T. For each wire u), 
with terminal qubits q{uj,0) and q{uj, 1), we define 

:= [q{uj,0) - qiuj,l)f , (7) 

=0 (1) for a satisfied (frustrated) wire. Then the error metric is 

w 

■= ^ ^UJ, (8) 

the number of frustrated wires. We do not engineer the error metric operator 
here. Solving the given problem requires minimizing the wire error subject to 
the triode relations. 

Each triode is associated with three orthogonal eigenmodes \9) defined so 
that the binary variable Sg has the value 1 for the mode |^) and in the other 
two modes. 

The three-dimensional Hilbert space of the triplet modes of triode r we 
designate by Ht- We define an auxiliary mode space of T disconnected triodes 
as the tensor product ^^^i 'Ht = W^, a Hilbert space of dimension 3"^. 

The total error form of the network is a lower bound on the number of 
wires that have to be changed to attain the solution. The number could be as 
great as W even if = 1 . 

We turn the error function into the Hamiltonian Hyj = ge^ diagonal in the 
qubit basis, with a coefficient g to provide the dimension of energy. 



5 



4 Computation model 



For qusa computation we take the effective Hamiltonian H^, + Hr{t), wliere 
Hr(t) is a small effective Hamiltonian term (possibly non-Hcrmitian) represent- 
ing relaxation processes that bring the network to the ground mode of Hw 
Hr{t) is discussed in the following. 

4.1 Continuous statistical projection 

To estimate qusa speed-up, we represent a development subject to statistics as 
a continuous projection of a dynamical development not subject to statistics. 

Consider a pair of identical protons 1 and 2. We freeze their spatial mode to 
a fixed antisymmetric wave-function tp-i2{xi, X2) so that only the spin degrees of 
freedom need be considered. The individual spin modes form two-dimensional 
Hilbert spaces Hi and H2- We consider: 

• an unsymmetrized tensor-product ("comparison") Hilbert space H' = 
Hi ® H2 ^ Hi (B Hs, the direct sum of the physical triplet subspace 

H = Ht and the singlet subspace Hs, 

• a particle exchange operator X12 : H' — > H' , 

• a projection P12 := ^(1 + X12) on the physical triplet subspace H = HtC 



• a symmetrized Hamiltonian if 12 : Ht ^ Ht invariant under proton ex- 
change: Xi2-ffi2 = -ffi2-'^i2, so that X12 is a constant of the motion. 

• the extension H of the spin Hamiltonian operator H12 from the triplet 
subspace Ht C H' to the entire space H', with H = on the singlet space 
for convenience. 

The orthohydrogen modes are symmetric in the spin variables. They vary 
over the effective three-dimensional Hilbert space of triplet modes Ht CH'. 

In the Hilbert space H', the time development of the orthohydrogen spins is 
then governed by a Hamiltonian H that maps the triplet subspace into itself: 



where |0) is a symmetric mode in H' . We now develop an equivalent represen- 
tation of the time development. The aim is to free the Hamiltonian from the 
mathematical conditions representing proton indistinguishability. 

We start from a symmetric initial mode |0) of H' at time and let it evolve 
for an infinitesimal amount of time dt according to a different comparison Hamil- 
tonian H': 



H' 



\dt) := (1 - iHdt) |0) , 



(9) 



\dt)' := (1 - iH'dt) |0) . 



(10) 
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H' is not subject to exchange symmetry, but we choose it so that its sym- 
metrization yields the actual Hamiltonian: 

Pl2H'Pi2 = H. (11) 

In general Idt)' is not symmetric under proton exchange. We restore particle 
indistinguishability by projecting Idt)' on Tit. This means symmetrizing Idt)' to 
form 

Pi2 \dt)' = (1 - iHdt) |0) . (12) 

The projection of the infinitesimal development (10) on Ht yields the actual 
development (9), up to higher order infinitesimals. We "continuously project" 
on Ht the development governed by H' . That is, we project after each interval 
At and take the limit At — > 0. This recovers the actual development generated 
by the symmetrized Hamiltonian H . 

This mathematical artifice of asymmetric time-development accompanied 
with continuous symmetrization permits us to estimate the speed-up due to 
quantum statistical computation. 

4.2 The speed-up due to statistics 

For qusa computation, we apply the continuous symmetrization of Section 4.1 
to the relaxation of the triode network; see for example (2). For a comparison 
network we work in the unsymmetrized tensor product Hilbert space Ti.' D 
7i^, suspending proton indistinguishability and removing all the consequent 
statistical relations. 

This means dropping the triode condition Qx + Qy + Qz — 1 for the weaker 
condition Qx + Qy + qz = ^ mod 2. The latter holds independently of statistics 
due to the composition of angular momentum alone, so it survives. 

usually has traps (local minima relative to all immediately adjacent 
energy levels) that slow classical annealing computation. We eliminate these for 
the comparison computation by redefining in a way that does not change 
the ground mode: 

i^u. = ge^ + g'e^ 51 { [1 - '?(^' 1)1' + [1 - '?(^' 2)]' + [1 - 9 3)]'} , (13) 

r 

where g(r, 1), ... are the three qubits of triode r. 

Since for the triode network each triode has only one qubit equal to 1, we 

have merely multiplied the previous Hamiltonian by ^1 + 2T^^ . This docs not 
change the ground mode. 

In the case of the EQU network, all the three qubits of a triode can be 1. 
If each frustrated network mode (e^ > 1) has a gradient toward the 

solution where all qubits are 1. This ground mode is quickly reachable by the 
EQU network even in classical annealing. 
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We define the effective Hamiltonians of the actual and comparison networks, 



H = H^+Hr, (14) 
H' = H'^+K. (15) 

Each describes a network with symmetric wire Hamiltonian Hyj or = Hy,, 
coupled to a heat bath by an effective coupling Hamiltonian respectively Hr 
and H^. H'^ is not symmetrized; its symmetrization yields H^. H describes 
the relaxation process of the actual triode network. H' describes the relaxation 
process of a comparison network obtained by replacing all pairs of identical 
indistinguishable spin 1/2 particles with pairs of distinguishable (e.g. non- 
identical) spin 1/2 particles. Correspondingly, all triode relations are replaced 
by EQU gate relations. 

We model the actual heat bath coupling Hr by coupling each proton spin u 
to a small Gaussian random time-varying magnetic field B at the site of that 
spin. B might be polarized along the principal direction x + y + z. We index 
the sites with the triode index r = 1, . . . , T and a binary index /? = 1, 2. While 
T enumerates the triodes (proton pairs), /? distinguishes the two protons in each 
triode. 

To preserve statistics we must demand that the two protons /3 = 1, 2 of each 
triode r experience the same magnetic field B{t). We may then write the actual 
relaxing coupling as 

Hr=gY,B{T)-g{T,p) (16) 

The comparison heat bath is a random magnetic field at each proton site. 
Unlike the actual heat bath coupling, the comparison heat bath coupling iJ'. is 
not invariant under proton exchange. Different protons in the same triode see 
different magnetic fields B{t, (3, ): 

H'^=gJ2B{T,P)-a{T,P) (17) 

r,0 

B{t, P) too might be polarized along the principal direction x + y + z. 

Let P be the symmetrization operator for all the relevant proton permuta- 
tions; it is not necessary to permute protons between triodes. We may arrange 
that the projected heat-bath coupling is the actual coupling, 

PH'^P = Hr, (18) 

by identifying the random magnetic field B(t) of the actual heat bath with the 
average of the two random magnetic fields of the comparison heat bath; 

BM.^fciii^ (19, 

The sum of two Gaussian random variables is also a Gaussian random variable. 
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Summing up, we have PH'^P = H^, PH'^P = Hr, and thus PH'P = H. 



Therefore we can apply the method of Section 4.1. 

The fact that H'^ — is already symmetric with respect to proton ex- 
change (as all its qubits are) does not introduce any unwanted constraint in 
the comparison relaxation process. It does not prevent the generation of tri- 
ode violations by H'^.. Hyj symmetry only reflects the composition of angular 
momentum, not statistics. 

Here we do not tackle the difficult problem of estimating the relaxation time 
of the actual triode network, or the comparison EQU network. We just know 
that H' will eventually drive the EQU network to its ground mode, and compare 
the relaxation times of the two networks. 

For this comparison, decompose the actual computation time AT into N = 
^ consecutive time slices Ai^ : ti < t < U+i of equal length At, with i = 
1,2, ...,N, ti = iAt. Take the relaxation within each Ati to be that of the 
comparison EQU network. At the end of each Ati project the network mode on 
the Hilbert space T-f" . Then take the limit At — > 0. This reproduces the actual 
computation. 

Let |0) be an initial symmetrical preparation, where all triodes must be 
satisfied and wires can be frustrated. Let this evolve into \t) at time t, with 
random phases as required. To consider the development of \t) inside the interval 
Ati, we resolve \t) as follows: 

\t) = \S,t) + \F,t) + \V,t). (20) 

\S,t) denotes a superposition of tensor products of triode eigenmodes (terms) 
which are solutions of the triode network (each with a random phase to rep- 
resent incoherence as necessary); its terms have satisfied triodes and satisfied 
wires. Most probably, relaxation randomly generates a solution with probability 
Ps{t) ■= {S,t\S,t) of the order of 1/2^ in po\y{Q) time. We assume this is the 
case at time ti. Since 

H^\S,t)^0, (21) 

we can assume that {S,t\S,t) remains approximately constant inside Ati. 

\F, t) is the component of \t) whose terms have satisfied triodes and at least 
one frustrated wire. Its probability pF(i) := {F,t\F,t) is initially close to 1 (see 
further below). 

\V,t) is the component whose terms have at least one violated triode; wires 
can be either satisfied or frustrated. Its probability is pv{t) ■= (^j^l^i^)- 

\V,t) is generated by the relaxation of the EQU network inside each Ati; it 
goes to zero with At and is annihilated by the projection at the end of Ati. 

We compare the rate of relaxation of the triode network to that of the EQU 
network as follows. 

Inside Ati, the evolution is that of the comparison network where all triodes 
are replaced by EQU gates. Therefore prit) goes down at the rate of relaxation 
of the EQU network, building up the mode \V,t). 

We are particularly interested in the "take off" of the solution probability 
Ps{t) from O (l/2'5) to 0(1), say to ps{t) = 1/10. During take off, we can 
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assume a constant rate of relaxation k of the EQU network with prit) ~ 1- 
That is, PF{t) changes from pFiU) ~ 1 at the beginning of Ati to 

PF{ti+i) « (1 - kAt)pF{ti) « (1 - kAt) 

at the end of Ati. Correspondingly, pv{t) changes from to about kAt. 

The projection at ti + Ati therefore reduces {t\t) by about kAt. Renormal- 
izing then multiphes ps{t) by about (1 — kAt)~^ f« 1 + kAt at each Ati. After 
a time NAt = AT we have 

ps{ti + AT) = psiti) lim^ (1 + kAt)^ = ps{ti)e'''^^. (22) 

Thus ps{t) becomes 0(1) in a time AT such that e''^^ = 0(2^). This 
means kAT = 0{Q), the number of qubits. The relaxation time AT of the 
actual network is comparable with that of the easy EQU network. 

Let us check that assuming a constant (average) relaxation rate does not in- 
troduce unwanted restrictions. We show that the continuous projection method 
works in the same way in presence of fluctuations of the expected energy of the 
EQU network. Assume this energy goes up in the time interval [ti,tj] , with 
tj > ti, namely {tj \ \tj) > {ti \ \ti) . On the basis of (22), this brings ps{t) 
down but, as readily checked, what is lost is exactly regained when the expected 
energy goes back to {ti \ \ti) at some time th > tj. 

Summing up: (i) We have defined an effective symmetrized Hamiltonian 
H for the relaxation of the "hard" triode network, (ii) We have defined a 
comparison non-symmetrical Hamiltonian H' whose symmetrization yields H: 
PH'P = H. (iii) Furthermore H' describes the relaxation process of an easy 
EQU network, obtained by replacing all triode relations by EQU relations, (iv) 
The relaxation of the hard triode network can be obtained by continuously 
projecting that of the easy EQU network on the symmetrical subspace H^. By 
eq. (22), this yields comparable relaxation times for the hard and easy networks. 

Satisfying relations implemented by particle statistics does not take time in 
annealing computation. 

If quantum annealing c:an rc^acli the solution 111... of a general EQU network 
in polynomial time, then in qusa computation NP=P. 

Qusa computation survives decoherence as well as general annealing compu- 
tation does. They both avoid this basic difficulty of reversible quantum compu- 
tation. 

5 Discussion 

5.1 The origin of the qusa speed-up 

What makes quantum computation more efficient than classical computation 
has been called a "most pressing" question for the advancement of the field 
(Mahler 2001). The speed-up of the quantum algorithms (Deutsch 1985, Shor 
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1994, Grovcr 1996, among others) stems from the fact that a quantum transi- 
tion is jointly influenced by an initial and a final extra-dynamical selection. It 
may be regarded as extra-dynamical in origin (Castagnoli & Finkelstein 2001). 
This kind of process is richer than classical computation, which is a dynamical 
development of one selection alone. 

Qusa computation can be seen as a dynamical development generated by a 
symmetric time- varying Hamiltonian H with a time-dependent magnetic field B 
(varying randomly for annealing). The extra-dynamical origin of its speed-up is 
H symmetry itself. This originates extra-dynamically in proton indistinguisha- 
bility, not dynamically. By the way, one can see from equations (10) and (12) 
that symmetrizing the network mode and symmetrizing H' are equivalent. 

It is as if proton indistinguishability provides for free an extra-dynamical 
symmctrization engine continously acting on an un-symmctrized H' or the de- 
velopment generated by it (relating the easy EQU network). The comparison 
with the quantum algorithms is easier if we consider the development. 

Qusa speed-up comes from the fact that the Boolean statistical relations are 
satisfied through the projection of a development that is unaffected by them. 
Correspondingly, we have seen that in each Atj there is a quantum transition 
jointly influenced by an initial and a final extra-dynamical selection. This is 
as in the quantum algorithms, but for the fact that qusa computation always 
selects a predetermined subspace, the one that satisfies statistics. 

5.2 Conclusions 

Qusa speeds quantum ground mode computation up by implementing the gates 
(or partial gates) of a Boolean network with statistical symmetries, and only 
the wires through an energy function which is zero when the wire is satisfied. 
Logical relations associated with statistical symmetries do not slow down the 
development of the initial mode toward a ground mode where both gates and 
wires are satisfied. This is unlike classical ground mode computation, where 
logical relations reduce relaxation rate. 

Qusa computation develops quantum parallelism through the incoherent su- 
perposition of parallel computation paths (i.e. through mixtures). It abandons 
the delicate superposition of coherent parallel computation paths in reversible 
quantum computation but introduces the almost indestructible superposition of 
different permutations of identical particles subject to a given statistics. This 
greatly reduces the problem of decohcrence. 

This synthesis of ground mode computation and quantum statistics appears 
to be a promising architecture for robust quantum computing. 

A natural next step in this research is to design physical systems exhibiting 
the qusa effect. 
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